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Two predictions are made for properties of the ferromagnetic superconductors discovered recently. 
The first one is that spin-triplet, p-wave pairing in sucli materials will give the magnons a mass 
inversely proportional to the square of the magnetization. The second one is based on a specific 
mechanism for p-wave pairing, and predicts that the observed broad anomaly in the specific heat 
of URhGe will be resolved into a split transition with increasing sample quality. These predictions 
will help discriminate between different possible mechanisms for ferromagnetic superconductivity. 
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Recently there has been a considerable amount of work 
on what is believed to be true coexistence between fer- 
romagnetism and superconductivity. In UGe2 0, ZrZn2 
0, and URhGe J.], superconductivity has been observed 
at very low temperatures inside the ferromagnetic phase. 

There have been several theoretical suggestions for the 
cause of this exotic superconductivity. Early theories 
proposed magnetically mediated spin-triplet supercon- 
ductivity in ZrZn2 Q, and recent band-structure calcu- 
lations have concluded that this is a possibility in these 
materials 0], and Machida and Ohmi have give 
arguments for the pairing to be of p-wave, non-unitary 
spin-triplet nature. On the other hand, a conventional 
phonon mechanism has been proposed in Ref. |^ and 
Ref. 13 has suggested that the superconductivity is due 
to coupled CDW and SDW fluctuations. It is not obvi- 
ous from these considerations why the superconductivity 
is observed only in the ferromagnetic phase. To explain 
this, Sandeman et al. have proposed a density-of-states 
effect, that exists only in the ferromagnetic phase, as the 
source of the superconductivity in UGe2 ^3 ■ Kirkpatrick 
et. al. have proposed an explanation for the ob- 

served phase diagram that is based on an enhancement of 
the longitudinal spin susceptibility in the ferromagnetic 
phase by magnons, or magnetic Goldstone modes. 

In this Letter we predict two observable effects that al- 
low to partially discriminate between these various the- 
oretical mechanisms. First, we show that if the order 
parameter (OP) is of non-unitary spin-triplet type, the 
ferromagnetic magnons, or transverse spin fluctuations, 
become effectively massive in the superconducting phase, 
with a mass that is inversely proportional to the square of 
the magnetization. This prediction depends only on the 
symmetry of the OP, and is independent of the mech- 
anism for superconductivity. The second prediction is 
based on our previously suggested mechanism for super- 
conductivity in these systems ^3 ■ It predicts that 



the specific heat anomaly observed in URhGe around the 
superconducting transition temperature actually con- 
sist of two distinct features that will be resolved as the 
sample quality increases. These two features correspond 
to the pairing of electron spins parallel and antiparallel, 
respectively, to the direction of the magnetization. In 
what follows, we will first give our detailed results and 
simple intuitive explanations, and then we will sketch 
their technical derivation 1^. 

General symmetry principles |l4| show that, in the 
presence of a nonvanishing magnetization in z-direction, 
the magnetic susceptibility tensor has the structure 
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with xt,± and xl the transverse and longitudinal mag- 
netic susceptibilities, respectively. In the absence of su- 
perconductivity, YT is massless, since it is the magnetic 
Goldstone mode |l4| . Spin-triplet superconductivity ex- 
plicitly breaks the spin rotation symmetry and leads to a 
mass in xt E3 • ^'^^ small frequencies fl and wavevectors 
k we find 
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Here Ct — fl/Aep and k — k/2k-p are the frequency and 
wavevector made dimensionless by means of the Fermi 
energy ep and the Fermi wavenumber k■p^ respectively. 
a+ and b are real constants on the order of unity, and a_ 
is imaginary with a modulus on the order of unity. 5 — 
(5/4eF, with 5 the Stoner splitting of the Fermi surface, 
which is proportional to the magnetization m. In a phase 
where only up-spin electrons are paired, with A| the up- 
spin superconducting energy gap, we find for the mass 
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FIG. 1: Circles denote the experimental i^Uc heat co- 
efficient of URhGe |3- The sohd hnes show the theoretical 
expectation if only half of the electrons pair (a), and if there 
are two transitions within 5% of one another (b). See the text 
for additional information. 
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Asymptotic analysis yields f{x — > oo) = Inx, and f{x 



0) = 0.85. 
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Such a mass in the transverse magnetic 



susceptibility is in principle observable by neutron scat- 
tering. An estimate using Af « 2Tc| = 1 K, 5 « eF/2 = 

10^ K, and fcp = yields 2kF^i = 3 x IQ-^ . 

This is roughly a factor of 30 smaller than the current 
lower wavenumber limit for inelastic neutron scattering 
Since n depends roughly linearly on A|, a higher 
value of Tc would help push fi into an observable regime. 

Our second prediction concerns the detailed structure 
of the specific heat anomaly oberserved in URhGe 
The experiment shows a broad peak with a maximum 
close to the resistive Tc, which is evidence for true bulk 
superconductivity. However, for low temperatures the 
specific heat coefficient does not go to zero, but rather 
approaches a finite value that is roughly half its value in 
the normal phase, see Fig. ^ Ref. y has interpreted this 
behavior as an indication for the pairing being restricted 
to one spin projection, while the electrons with the op- 
posite spin projection retain a Fermi surface down to the 
lowest temperatures measured. With a relative specific 
heat discontinuity AC/C adjusted to fit the data, and a 
temperature dependence below Tc given by C{T)/T oc T 
|l7j . this suggests that, with improving sample quality, 
the data should approach the solid curve shown in Fig. 
^a). This is not the only possible interpretation. We 
have extended our previous calculation of the supercon- 
ducting Tc to include both up-spin and down-spin pair- 
ing. For the superconducting gap, which is a 2 x 2 matrix 
A in spin space, we assume )18l] 
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with k a unit wavevector. This is the OP of the A2 
phase in He^ |l9| . For generic parameter values we find 
that there is only a small difference between the values 
of and T^i in the ferromagnetic phase, on the order 
of 5-10%. This is comparable to the difference Fay and 
Appel found within a simpler theory that could not 
describe the strong enhancement of Tc in the ferromag- 
netic phase. Consequently, one expects that in a perfect 
sample the specific heat data would qualitative look as 
shown in Fig.^b). In this interpretation, the broad peak 
hides a split transition, and the low-temperature tail of 
the specific heat coefficient is due to non-superconducting 
parts of the sample. (See, however, Ref. |2(J.) It is inter- 
esting that such a split transition has been observed in 
the heavy-fermion superconductor UPta [21I (which 
is not ferromagnetic, and where the splitting has a dif- 
ferent physical origin), but only after a long period of 
improving sample qualities. 

We now turn to the origin of these results, starting 
with the mass in the magnon. Any spin-triplet OP can be 
represented as a vector 4> in spin space that is isomorphic 
to the 2x2 spin matrix representation For the OP 
given in Eq. Q), this vector is proportional to 

0== (At - A^,i(At + A^),0) . 

We see that </>* is not parallel to 4>- This an example of 
a "non-unitary" OP, and Ref. [3 has argued that any OP 
in ferromagnetic superconductors is likely of this general 
type. In a Landau theory, the lowest order term in the 
free energy density that couples the magnetization vec- 
tor M = (0, 0, m) and the superconducting OP that is 
allowed by symmetry is thus of the form 
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with c a real constant. This shows that the presence 
of a nonzero superconducting OP, A| is equivalent to 
an effective magnetic field hes ot. A|. General symme- 
try considerations [2^ show that the transverse mag- 
netic susceptibility must then have a mass that is given 
by Xr^i^T^) — hcs/m. This explains the existence of a 
mass, and the dependence of /i^, Eq. (p?ajl . on A|. 

For an understanding of the remaining structure of 
/i^, as well as of our second prediction, we return to 
our previous work on the superconducting transition of 
the up-spin electrons in that paper, we con- 

sider a model of free electrons with a static, point-like 
spin-triplet interaction with amplitude Ft. For super- 
conducting OP fields we choose, !Fi{x,y) = tjji {x)ip'^ (y) 
and J-i{x, y) = ip \,(x)'ip lijj) , with 'ipa-i^) an electronic field 
with spin index a and space-time index x. The OPs, i.e., 
the expectation values {J-^{x,y)) — F^{x — y), are the 
anomalous Green functions. 

Using a technique similar to that employed in Ref. 
we have derived coupled equations of motion for the Fa-, 
the normal Green functions G^, and the magnetization 
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m, in a systematic loop expansion. Because our pairing 
mechanism is due to magnetic fluctuations, one needs to 
go to one-loop order in order to obtain superconductiv- 
ity. To that order, and for p-wave orbital symmetry, we 
find the following Eliashberg-type equations for the gap 
functions A^- and the normal self-energies So-, 

A.(fc) = r* f XL{k-q)^M)/da{q) , (6a) 

Jq 
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with a ='f , 1= -(-,—. Here we use a wavevector- 
frequency four- vector notation k = {)<i,iuJn)^ and 
Sq = ^E„(V^)Eq- We have introduced d„{k) = 
G-\k) G-\-k) + ^, (A:) At (fc), and the G~\k) = ic^„ - 
(fc^/2me — ^l — cr5) — Scr(fc) are the inverse 'normal' Green 
functions, where 5 = Ttm is the Stoner energy. The mag- 
netic equation of state to one-loop order is given by, 

m = ^T^r{j3g) + Y / dxdy^Xjkix -y) 

xtr [73e(-a;)7je(a; - y)7fca(2/)] • (7) 

It is easy to show that the magnetic susceptibility x 
that serves as the effective potential in these equations is 
proportional to the physical spin susceptibility, Eq. 
In a Gaussian approximation j2J| we find 

p 

Xij^{x-y) = 5,j S{x~y)+Y tr [Giv ' x) H Q{x - y) ] . 

(8) 

In these equations, the (70,7) = (ca ctq, (t^ ® (Ji,aQ ® 
CT2 , fa (8) (Ta ) are generalized Nambu matrices, and the Q 
are 4 x 4-matrix Green functions that represent the nor- 
mal and anomalous Green functions in this basis. Tr 
denotes a trace over all degrees of freedom, and tr de- 
notes a trace over the discrete degrees of freedom. 

An evaluation of Eq. (jS]) leads to a transverse suscep- 
tibility tensor that has the form of Eq. with a^. = 1, 
a_ = — i, b = 2/3, and jj? as given in Eq. l|Sa|l . The gen- 
eral Landau theory argument given above, see Eq. |(3J), 
shows that this form of x is generic, and not an artifact 
of our Gaussian approximation. (See, however, Ref. B) 
The explicit calculation also shows that the prefactor of 
the A^ in the mass is a singular function of the magne- 
tization and the temperature. This implies that a strict 
Landau expansion in powers of both the magnetic and the 
superconducting OPs is singular, since the coefficients of 
the Landau function do not exist. This breakdown of the 
Landau expansion is due to the soft particle-hole excita- 
tions in an itinerant electron system, and is analogous to 
an effect that has been discussed for quantum ferromag- 
nets in the absence of superconductivity (25| . 



Our second prediction hinges on a solution of the gap 
equations, Eqs. ©, which require the magnetic suscep- 
tibility tensor x as an input. A complete numerical solu- 
tion of these strong-coupling equations would very diffi- 
cult, and does not seem worthwhile in the absence of de- 
tailed experimental information about x- I^^f. ^3 we 
emphasized that the relative values of Tc_| in the para- 
magnetic and ferromagnetic phases, respectively, can be 
obtained within a simple McMillan-type approximation, 
even though the absolute values obtained by this method 
are probably not reliable. Here we adopt the same phi- 
losophy, with the aim being to compute the relative mag- 
nitudes of both Tcj and T^i- In the paramagnetic phase 
this is straightforward, using the Gaussian approxima- 
tion, Eq. JHl, for the magnetic susceptibility. We find 
results for Tcj and Tc| that are comparable to those ob- 
tained by Fay and Appel Q, i.e., a ratio Td/Tci that is, 
for generic parameter values, on the order of 1.1. 

In the ferromagnetic phase we need to take into ac- 
count the effect discussed in detail in Ref. 0, namely 
an enhancement of the longitudinal susceptibility due to 
a coupling of the latter to the ferromagnetic Goldstone 
modes |l26|. This requires a calculation of x to one-loop 
order, beyond the Gaussian approximation given in Eq. 
lIHll . For Tc| this calculation has been performed in Ref. 
I12L For Tci there is the additional complication that 
inside the superconducting phase, x depends on the up- 
spin gap A|. This introduces a feedback effect that is 
characteristic of any purely electronic mechanism for su- 
perconductivity. Since the overall effect of this feedback 
is rather small (see below), we can approximate the tem- 
perature dependent A| in x by 2Tc| . We have found that 
Tcj, like its Tc| counterpart, is enhanced over its value 
in the paramagnetic phase by a factor of 50-100, and for 
generic parameter values it is within about 10% of TcT- 
For the specific heat this leads to the prediction shown 
qualitatively in Fig. ^ namely, two transitions in close 
proximity. Our calculations therefore suggest that the 
specific heat peak measured in URhGe actually consist 
of two unresolved peaks. We note that the width of the 
experimentally observed peak easily encompasses both 
mean- field discontinuities. In contrast to our first predic- 
tion, the second one is specific to the pairing mechanism 
we have considered. If experiments on cleaner samples 
should show no indications of two closely spaced transi- 
tions, that would be a strong argument against the spin- 
fluctuation induced pairing we have considered. 

We finally discuss our second result in some more de- 
tail. The salient point of our mechanism for a strongly 
enhanced superconducing in the ferromagnetic phase 
is the coupling of the massless (on the superconducting 
boundary, where A^ = 0) transverse spin susceptibility 
to XL via mode coupling effects. Because xl serves as 
the pairing potential for the A| ordering, this suggest 
that in the superconducting phase the pairing potential 
decreases, since xt decreases for all wavenumber and fre- 
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quencies with increasing A^. For spin-triplet supercon- 
ductivity in a magnetic field one generically expects an 
ordering of the spin down electrons, A| ^ 0, at some 
temperature Tc| < T^^. Since the pairing potential for 
A| is also XL, this implies there is a weaker tendency for 
the onset of A| given A| 7^ 0, than there is for the onset 
of A| ordering. Another general mechanism decreasing 
Tcj compared to Tc| is that the density of states at the 
Fermi surface for the down-spin electrons decreases with 
increasing magnetization. Naively, one therefore expects 
that for small magnetization the cutoff in Eq. ||2Jl, which 
is proportional to greatly suppresses the A| order- 
ing, while for large 8 the density of state effect should 
decrease it. This would result in, for example, a spe- 
cific heat discontinuity at Tcj , and a substantial residual 
specific heat coefficient at lower T, since the down-spin 
electron would not undergo a superconducting transition 
until perhaps immeasurably low temperatures. This is 
certainly consistent with the experimental results in the 
URhGe system, and it is the interpretation given in Ref. 
I3- However, our detailed numerical work suggests that 
this is not the generic situation for the pairing mecha- 
nism we have considered. Clearly, with improving sam- 
ple quality that will sharpen the feature in the specific 
heat, this question can be decided experimentally. 
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